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Abstract. Using an asymptotic phase representation of the particle density operator 
p{z) in the one-dimensional harmonic trap, the part 5pp(z) which describes the Friedel 
oscillations is extracted. The expectation value (<5/5p(z)) with respect to the interacting 
ground state requires the calculation of the mean square average of a properly defined 
phase operator. This calculation is performed analytically for the Tomonaga-Luttinger 
model with harmonic confinement. It is found that the envelope of the Friedel 
oscillations at zero temperature decays with the boundary exponent v = (K+l)/2 away 
from the classical boundaries. This value differs from that known for open boundary 
conditions or strong pinning impurities. The soft boundary in the present case thus 
modifies the decay of Friedel oscillations. The case of two components is also discussed. 
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1. Introduction 

Recent experimental successes in obtaining degeneracy in three-dimensional ultracold 
Fermi vapors [1-7], possibly in combination with microtrap technology [8-13], make it 
conceivable to realize the quasi one-dimensional neutral Fermi gas confined in a trapping 
potential without the complications due to contacts and impurities. 

Friedel oscillations are a principal feature of a degenerate Fermi gas when 
translational invariance is broken. Usually, impurities [14] are the cause. However, 
boundaries can also be responsible for Friedel oscillations. The spatial period of 
Friedel oscillations is = ir/kp, where kp is the Fermi wave number. This effect 
is particularly pronounced in one spatial dimension because then the susceptibility 
becomes logarithmically singular at 2kp due to perfect nesting. 

It is known from the theory of one-dimensional fermions confined between hard 
walls (bounded Luttinger liquids = BLL) that interactions modify the decay of Friedel 
oscillation away from the boundary [15-17]. 

In this article, we investigate the Friedel oscillations at zero temperature both for the 
spin-polarized one component system, when s-wave scattering is forbidden, as well as for 
the two component system. We apply a recent model of N ^> 1 interacting fermionic 
atoms in one spatial dimension confined in a harmonic trap [18]. The model which 
can be termed "Tomonaga-Luttinger model with harmonic confinement" is analytically 
solvable in a similar way as the Luttinger model (cf. e.g. [19-22]) using bosonization. 
The present model is simpler than the Luttinger model in that it has only one (non- 
chiral) branch in contrast to the two chiral branches of the latter which result from an 
artificial split of an otherwise continuous band. 

It is found that near the classical boundaries, Friedel oscillations decay in a way 
which differs from the known result for BLL [15-17]. It also differs from the result for 
an infinitely strong pinning impurity [23] which acts as an invariant hard wall under 
scaling [24]. 

The calculations become possible because the fermion density in the harmonic trap 
can be decomposed asymptotically (i.e., for N ^> 1) into a slowly varying part and a 
part describing the Friedel oscillations. Both parts involve a specific phase operator for 
which a free field theory is available. 

Friedel oscillations in realistic Fermi gas are difficult to observe experimentally 
at least for two reasons: the integrated total mass in the one- dimensional Friedel- 
oscillations is of the order of one atom (cf. e.g. [25]) though repulsive interactions 
increase their weight. Furthermore, temperature effects blur the oscillations unless 
ksT < hue where u>e is the longitudinal trap frequency [26]. Provided these exceedingly 
low temperatures can be realized, we can adopt the argument in [25]: it is conceivable 
to use an array of short micro traps each filled with a reduced number of atoms (thus 
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avoiding instabilities). The oscillations within each trap then add up and lead to a total 
effect that is enhanced by the number of traps. Using micro fabrication techniques it 
should be possible to combine 100 traps on one substrate leading to a signal that may 
become within reach of advanced imaging techniques. 



2. Theoretical framework 

2.1. Tomonaga-Luttinger model with harmonic confinement 

Here, we give a short review of the model which is used for the present calculation. 
More details are given in [18]. 

We start from an effective pair interaction 

V = \ E V{m,p;q,n) {c+c q )(c+c n ) (1) 

mnpq 

between one-dimensional spin polarized fermions. The fermionic creation and 
destruction operators c+ and c q are taken in the basis of harmonic oscillator wave 
functions. Thus the harmonic trap topology is exactly represented. 

The Hamiltonian for interacting fermions considered here follows uniquely by 
retaining those parts in the fermionic pair interaction operator which are expressible 
in terms of density fluctuation operators p(m) = J2 P c p+ m £p- These are 



V(m,p; q, n) -> V a (\q - m\) 5 m _ g , n _ p + V b (\q - m\) 5 ? _ m , n _p + V c (\q - p\) 5 m+q>n+p . (2) 

The interaction matrix elements V a (m), V&(m), and V c (m) correspond to the 
Luttinger model coupling functions g^ip), gi(p), and gi(p), respectively. V a and Vf, 
describe forward scattering and V c describes 2k f (backward) scattering. In [27], it is 
shown that the retained matrix elements are dominant in the limit of large N. This is 
related to approximate momentum conservation during collisions in the trap. 

In the next step, the linear dispersion of free harmonic oscillator states and 
the addition of the anomalous vacuum is utilized to bosonize the original fermionic 
Hamiltonian in terms of canonically conjugate Bose operators d and d + in accordance 
with Kronig's identity [28]. This gives the bosonic form 



H = \h^Y. m {dJ m + d + J m ) -\Y. V cM m (d 2 m + di 2 ) (3) 

Z m>0 Z m>0 

+ \ E V d m ) [tkm + dtm] 

/ m>0 

+ \ ^ V a (m) m (dj+ + d + J m ) + \ E V b (m) m (d 2 m + rf+ 2 ) . 

Z m>0 Z m>0 
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of the Hamiltonian. 

Equation (3) is the version for a one-component Fermi gas, i.e., for spin polarized 
fermions. Usually, the remaining "p-wave" interaction is small. However, it has been 
demonstrated recently that Feshbach resonances can make it relevant [29]. 

Our simplified interaction Hamiltonian is integrable. In the sense of the Luttinger 
liquid phenomenology, we expect that the boundary exponent for Friedel oscillations at 
zero temperature and for N ^> 1 is invariant against details of the interaction. It is also 
stressed that the trap potential is exactly incorporated. 

In the one-component system, backscattering dominates as demonstrated below. 
Accepting this, the validity of our approach could be verified analytically by perturbation 
theory in the fermionic Hilbert space as well as by exact numerical diagonalization of 
the fermionic problem [18]. 

2.2. Dominance of backward scattering in the one- component gas 

The existence of only one branch results in restrictions on the values of the interaction 
coefficients V a and V& in the one-component case: their contribution is small in 
comparison to V c [27]. This can be demonstrated analytically by using WKB wave 
functions (cf. equation (21) below) in the calculation of V(m,p;q, n): starting from an 
effective "p-wave" potential in one dimension 



V(z)=V p a 3 p d 2 z S(z), 



(4) 



one obtains for (m, n, p, q) 



O(N) > 1 



V(m,p;q,n) = / dzdz' i\) m {z)i\) q {z) V(z - z') i/j p (z')ijj n (z') 



(5) 



f 4 

= V p al I dzi) m (z)i) q (z) d\ {tp p (z)ip n (z)} -> - 




F(s), 



with 



F(s) 



COS 2 (7Ts/2) 



S 



III ■ (j j) II. 



(6) 



Thus, each individual backscattering term in (2) belongs to the dominating s = 
contribution while only some terms of V a ^-type do this. Thus backscattering dominates. 
The interaction coefficient V c = V c (l) can be expressed as 



2. 3. One-particle operator and phase fields 

The third contribution on the r.h.s. of equation (3) represents a one-particle operator 
V\. It originates from rearranging operators in equation (1) to bring backscattering into 
a bilinear form of density fluctuation operators. The one-particle operator is neglected 
in the BLL as pointed out in [30]. It does not alter boundary exponents, but has 
quantitative effects on other properties. V\ is exactly taken into account in the present 
model. 

The central dimensionless coupling constants and the renormalized level spacings 
for the model, equation (3), are given by [18] 



hwt + V a (m) - (V b (m) - V c (m)) 

hue + V a {m) + (V b (m) - V c {m)) ' 1 j 

e m = y/(hu t + V a {m)y - (V b (m) - V c {m))\ 

respectively. 

For simplicity, the dependence of K and e on m is occasionally suppressed. Note 
that V a (m) — > V a implies that there are no interaction effects due to this matrix element 
in the one component theory and V a is strictly zero while V b can be neglected, i.e., the 
central coupling constant is 

with V c given by equation (7). 

In some physical quantities, especially in the mean square phase fluctuation 
calculated below, the neglect of the m-dependence leads to inconsistencies because 
stability requires the coupling constants K m to approach unity for large m [20]. In 
analogy to the Luttinger model, we write approximately 




K m = 1 + (K - 1) exp(-mr). (10) 
An estimate for r is R/ L F <^ 1 where R is the spatial range of the interaction. 
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Assuming the same exponential decay for the interaction matrix elements 



V c {m) = V c e~ mr , r«l, (11) 
the one-particle operator can be rewritten as 



l r 

v ^T^L dn 

with the phase field 



g— r+2m g— r— 2iu 



+ 



]_ g— r+2iu J g— 2iw 



d u 0odd(u), (12) 



odd ( M ) = - ((/>(u) + + H - <f>(-u) - <p + {-u)) = W sin(rm) (d n + d+) . (13) 

1 n=l » n 



and is the bosonization phase operator [31], 



CO 1 

( W )E-^-e^ +, ^4 (14) 



n=i 



In diagonalizing equation (3), the one-particle operator contributes a c-number shift 
which is taken care of by the new phase operator 



$(u) = 0«w(u) + &(u), ( 15 ) 

with 

««)=^(B^)- (16) 

In contrast to o dd, the new phase field <E> is homogeneous in the operators / and 
/ + which diagonalize the Hamiltonian: H = J2 m me mfmfm + zero mode contributions. 
The zero mode plays no role in the present context. The representation of the phase 
field $ in terms of the diagonalizing operators is 



CO 



*(«) = E V — SinmU (/™ + /m)- ( 17 ) 



m=l 



3. Decomposition of the density operator for the harmonic trap 

The density operator p{z) in Fock space for fermions in the harmonic trap is 
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m=0,n=0 



(18) 



Bosonizing p(z) using the auxiliary field method in [31], i.e., 



C m C n 



2lx r 2n dudv 



JO 



4vr 2 



,i(m«- nv) 



-i(N-l)(u-v) 
I — e —r)+i(u-v) 



(19) 



x exp j— % (J) + (u) — (f) + (v ))} exp |— % (J)(u) — 0(f)) | 



implies a projection onto the subspace of N fermions: p(z) — > Pn{z). 

We will extract from Pn(z) a slowly varying part and a part associated with the 
Friedel oscillations. 

After bosonization, Pn(z) is given by 



/7T 
-7T 



* dudv e -i(tf-i)(«-«) 



4-^2 ]_ _ e -r;+i(u— u) 



e -j^+(u)+i^+(ii) e -i<^(M)+i^(u) 



X 



^ i) m {z)i) n {z) e 



imu—inv 



ra=0,n=0 



(20) 



Here, 77 is a positive infinitesimal. 

Using the form factor Z(z) = yjl — z 2 / L 2 F of the harmonic trap, where 2L F = 
2£^2N - 1 is the semi-classical extension of the Fermi sea, the harmonic oscillator 
wave functions ifj m can be approximated for N ^> 1 by their WKB form (a = l/£, I: 
oscillator length) according to 



ipm(z) 



2a 2 



1/4 



ir 2 mZ 2 (z) 



cos 



dxk m (x) - . 



(21) 



Because of the rapidly oscillating phase factor exp(— i(N — l)(u — v)), an expansion 
around the Fermi edge with 



k m {z) = ay/ 2m + 1 — a 2 z 2 w k F Z(z) + 



m 



L F Z{z) 



, k F = aV2N-l, m = N-l + m{22) 
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and 



J dxk rn (x) ss \zZ{z) + L F arcsin j + m arcsin (23) 



kpz(z) + m arcsin 



z 



Xp, 

is reasonable because it will compensate that factor. Explicitly, z(z) is given by 



z z 
2k F z(z) = k F Z(z)z + k F L F arcsin — = k F Z(z)z + (2iV — 1) arcsin — . (24) 

Lip hp 

We consider the sum over m in equation (20). Extending the m-summation to — oo 
and setting m = N in phase insensitive terms, leads to the asymptotic expansion 



oo / Ott 2 ™ 2 \ 1/4 

E Ms) ^^\ms) SN ' 1)U (25) 

x | e ^5(,)-^(iv-i)/2 ^ _ A^J V4 6{u + Mq(2)) 

It is noted that the core states are not properly represented in the expansion (25). 
As a consequence, one does not obtain correct information about the operator for the 
average density. Thus we will retain only the fluctuating part — > Sp which is of main 
interest here. 

Applying this expansion to equation (20), we find to lowest order in d u /N 



Sp(z) = - d z (j) odd (u (z)) cos2[k F z(z) + 4> o dd(u (z))]. (26) 

71 7T1]LpZ{Z) 

The phase field 4> dd{ u ) is found to be identical to that of equation (13). This phase 
operator (f) dd plays a central role in the present investigation. 

The density operator (26) consists of two parts involving the phase operator (j) oc id, 
a gradient term for slow spatial density variations and a rapidly oscillating term 8p F 
describing the Friedel oscillations. This structure is well known from the theory of the 
Luttinger model. In the latter case, the argument of the phase operator is the spatial 
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coordinate. In the present case, however, a non-linear relation between the spatial 
coordinate z and the variable u (z) in the phase operator according to 



appears which reflects the harmonic trap topology. Furthermore, it is a priori not clear 
what the right structure of the phase operator in the confined case is. Our calculation 
gives the answer in form of equation (13). 

4. Calculation of Friedel oscillations and boundary exponent 

In order to calculate the Friedel oscillations (5pf(z)) in the interacting ground state, 
we apply the Wick theorem in the well known form 



utilizing the homogeneous structure of the phase operator $ equation (17). Thus 



(cos2[k F z(z) - b(u (z)) + $(u (z))])o = cos2[k F z{z) - b{u {z))}e- 2{ ^ uo{z))2)o . (29) 




(27) 




(28) 



The mean square average with respect to the ground state becomes 




(30) 



which leads to 




-(K - 1) lnr 




(31) 




Here, the abbreviation 




(32) 



is introduced. Finally, one obtains 
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cxpj 2{<S>{u) M-vr 2 [(1 + e _ 2v)/2 _ e _ v cos(2m)] i/ 2 • W 



Considering r 1, 77 — > 0+, and 



2 ^ 



cos2u (^ = 2^— J -1, (34) 
leads to the result for the Friedel oscillations in the limit L F — \z\ 3> r 2 L F /8 

. (-1)" /rN*" 1 co S 2M( g )-6(iio(^))] 

<<w * ))o = - m; I2J • (35) 

It is seen that attractive interactions, K > 1, decrease the Friedel oscillations, while 
repulsive interactions, K < 1, increase them at any fixed position \z\ < L F . 

The backscattering phase shift in equation (35) is mostly a small correction as is 
seen by comparing equations (16) and (24), except when the coupling is exceedingly 
strong. 

It is noted that the coupling constant K goes to unity and b(u (z)) to zero for 
vanishing interactions. This answers a question put forward in [25]: the divergence of 
the envelope near the boundaries for free Friedel oscillations in the harmonic trap has 
the conjectured boundary exponent Kq = 1. For \z\ <C L F , i.e., well inside the trap, 
one obtains the free Friedel oscillations 



(-1)^ 

(5pf(z))oo = ~ „ T cos2k F z, (36) 
ZnL F 

in accordance with a corresponding result in [25], which was obtained by an asymptotic 
expansion of the exact particle density of non-interacting fermions. 

Returning to the interacting case and specializing to the region near the classical 
boundaries, where 



Z\z)^ 2(1 -\z\/L F ), (37) 

the boundary exponent for the decay of Friedel oscillations away from these classical 
turning points is 

v = (K + l)/2. (38) 
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In contrast, the corresponding result for BLL in the case of one component is 
v bll = K [15-17]. The latter value is also obtained for an infinitely strong pinning 
impurity which acts as a hard wall [23], a case when pinning strength does not scale [24]. 
The present soft boundary thus causes a slower decay for repulsive interactions in 
comparison to the BLL while the situation is reversed for attractive interactions. 

4-1. Dependence on trap parameters 

We discuss the question how trap parameters influence the value of the boundary 



exponent and the amplitude of Friedel oscillations. In view of k F = \Jm.A0Jt{2N — l)/h 
and equation (7), the central coupling constant K according to equation (9) depends 
on longitudinal trap frequency u>e and particle number N (apart from interaction data). 
Two cases will be considered: 

i The trap is made shallower, i.e., u>e is decreased keeping N constant. In this 
case, interactions become irrelevant and K goes to unity. Because the amplitude 
in equation (35) is proportional to 1/L F oc ^/uTi the Friedel oscillations vanish 
everywhere inside the huge trap as expected. However, they still increase towards 
the boundaries because the trap topology persists for all non-zero u>£. 

ii A kind of thermodynamic limit [32] is defined by making out oc 1/N and N large 
such that the Fermi wave number and hence V c and K ^ 1 remain constant. 
Again, the prefactor in equation (35) suppresses the Friedel oscillations everywhere 
inside the trap. Interestingly, at large but finite N the amplitude of Friedel 
oscillations increases faster towards the boundary when interactions are attractive 
in comparison to the repulsive case. 

5. Boundary exponents for two components 

We develop the corresponding theory for two components, e.g., two different hyperfine 
components of the same trapped fermionic atoms. We assume equal masses and equal 
trapping frequencies. The latter assumption is an approximation in the hyperfine case. 
The local densities for mass and composition: 




CO 



n( z ) [Cm+ + ^m— ^ 



(39) 



m=0,n=0 




n— 



m=0,n=0 



are treated as described in Section III. Two odd phase operators corresponding to 4>oM 
appear. They are defined by 
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in terms of the two bosonization operators for the components a — 1 and a = —1 



(40) 



0» = -i E "F= e jm (" + ^ 2 ) 6 mCT . (41) 
m=i V™ 

The relation of the 6-operators to mass and composition fluctuation operators d is 



The final representation of the total density operator is 



(42) 



Sp(z) = -d z (4>+,odd(u (z)) + 4>- )Odd (u (z 



7T 



(43) 



(-1) 



A? 



COS 



TTT]L F Z(z) 

Defining the basic phase fields 



2k F z(z) + 2(f) +jOdd (u (z)) + cos 2k F z(z) + 2<j>_ j0dd (uo(z)) 



®u(u) = —j= {<j> +j0dd (u) + IS <j)_ >odd (u)} + K(u) 



(44) 



with v — 1 for mass and v — — 1 for composition fluctuations, brings equation (43) into 
the form 



5p(z) = —d z (®i(u (z)) - &i(«o(*))) 

71 



(45) 



2(-l) 



A? 



COS 



7TT]L F Z(z) 

The analogue of equation (16) is 



2fc F z(z) +\^($i(u ( 2 )) - 6i(«o(^)))l cos(«s/2$-i(«o(*)))- 



v 7 4ei \l- e - r " 2m 

The quantity 6_! vanishes identically. 



(46) 



Coupling constants and renormalized level spacings are given by 
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K = 



and 



vVaJjnj) - \Vb\\{m) + vVbAjn 
\J (Hue + V a \\(m) + vV a ±(m)) + [V h \\(m) + uV b ±(m) - V c \\(m)} ' 



(hu e + V a \\ (m) 



V c \\(m)] 



(47) 



(hue + V a{l (m) + vV a± (m)) 2 - [V H (m) + uV b± (m) - V 4 (m)f 



nl/2 



(48) 



respectively. The subscript || refers to interactions between fermions of the same 
component while _L stands for different components [18]. 
Similarly, the composition operator reads as 



V2 . 
5a(z) = — d z ^^(u (z)) 

71 

2(-l) N 
H — , - ; sin 



nr]L F (z 

The central relation equation (15) generalizes to 



(49) 



2k F ~z(z) + y/2($ 1 (u (z)) - 6i(uo(^)))l sin(v / 2 $-i(mo^)))- 



oo rr 

$ V (U) = E\/- e ~ nV/2 sin ( UU ) {dnu + d+ v ) + K(U) 
n=l » ^ 



(50) 



E 

n=l 



n 



sinnM(/ ni , + /+ jy ). 



In order to be able to evaluate the phase fluctuations in the way given above, i.e., 
for free phase fields, the interaction coefficient V c j_ for backward scattering between the 
two components must be zero. This is usually not the case. Relying on insight from the 
Luttinger model [21], one can expect that for V c \\ > \V C ±\ the coupling V c ± scales to zero 
at low energies and that for V c \\ = V c ± > backscattering becomes irrelevant. In that 
case, K*_ x = 1 holds. This corresponds to spin isotropy in the Luttinger model. However, 
the present one-branch model hardly allows such a case in reality due to suppression of 
s-wave scattering in the parallel channel. With this provision, we give the results for 
two components based on the assumption that the field is a (renormalized) free 

field: 



{Sp F (z)) = 



2(-l 



TirjL F Z(z) 



-($ 2 _ 1 ( Uo W))o 



cos 



{2k 



F zlz 



V2bii 



u {z 



-<*?(«o(*))>0 



(51) 
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Note that no Friedel oscillations exist in the composition part because (5a(z)} = 0. 
With the renormalized value of the coupling constant K_ x — > K*_ x , the Friedel 
oscillations are 



(Sp F (z)) 



-1) 



N 



■ COS 



7lL F 

giving a boundary exponent 



(52) 



(53) 

4 2 v ; 

Again, this is different from the corresponding exponent 



K + if* 
2 

of the BLL [15-17]. 

The result equation (53) also applies to spinfull one-dimensional fermions [y = 
1 — > c: charge degrees of freedom, z/ = — 1 — > s: spin degrees of freedom) and is thus 
applicable to harmonically trapped one-dimensional electrons. 

6. Summary 

We have calculated the quantum interference phenomenon of Friedel oscillations for the 
Tomonaga-Luttinger model of interacting one- dimensional fermionic atoms trapped in a 
harmonic potential. We used bosonization techniques and an asymptotic representation 
of the density operator in terms of a particular phase field. The result shows that 
the boundary exponent for the decay of Friedel oscillations away from the classical 
boundaries differs from the result for bounded Luttinger liquids where the fermions are 
confined between hard walls. 
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